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We compute the tunneling density of states of doped multi- 
wall nanotubes including disorder and electron-electron in- 
teractions. A non-conventional Coulomb blockade reflect- 
ing nonperturbative Altshuler-Aronov-Lee power-law zero- 
bias anomalies is found, in accordance with recent experimen- 
tal results. The presence of a boundary implies a universal 
doubling of the boundary exponent in the diffusive limit. 

PACS numbers: 71.10.-w, 71.20.Tx, 72.80.Rj 



Carbon nanotubes provide a remarkable and excit- 
ing arena for mesoscopic transport phenomena involv- 
ing strong electron correlations fl. Recent experiments 
on single-wall nanotubes (SWNTs) have established the 
ballistic nature of SWNT transport and revealed the Lut- 
tinger liquid (LL) behavior of ID interacting fermions 
However, the situation is less clear for multi-wall nan- 
otubes (MWNTs), which are composed of several (about 
ten) concentrically arranged graphite shells, where ex- 
periments seem more consistent with diffusive transport 
e.g. showing typical weak localization features in the 
magnetoconductance. 

Quite remarkably, several experimental observations 
for MWNTs do not seem to fit into the framework of 
existing theories. The most prominent example concerns 
the pronounced power-law zero bias anomaly of the tun- 
neling density of states (TDOS) found at low energy 
scales, JS < 0.1 eV ffl, 



'(E) ~ E a 
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Furthermore, the TDOS at the end of the MWNT, while 
still of power-law form, is characterized by the boundary 
exponent a on d ~ 0.6. In this Letter, we show that these 
results can be understood in terms of a particularly ef- 
fective and non-conventional Coulomb blockade (CB) for 
tunneling into a strongly interacting disordered metal. 
Charge propagation on the MWNT is then effectively 
two-dimensional (2D), but for very low energy scales not 
probed in Ref. |j, E < E* with E* in Eq. @ below, a 
crossover to ID pseudo-gap behavior |J is expected. 

The key ingredient in CB theory is the probability 
P(E) that a tunneling electron excites electromagnetic 
modes with energy E in the system 0. The theory is 
meaningful if these modes are harmonic, and then P{E) 
directly determines the TDOS, 

HE) _ f°° dE , 1 + exp(~E/k B T) 



where vq is the non-interacting DOS. The probability 
P(E) is the Fourier transform of exp[J(£)] , with the phase 
correlation function (h = 1) 
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for a spectral density I(ui) of electromagnetic modes. For 
simplicity, we now focus on the zero-temperature case. 
Provided I(u>) remains finite for low frequencies, Eqs. (2) 
and (||) then straightforwardly lead to the power law (1) 
with exponent a = I(u> — > 0). 

In a conventional CB system, I(lo) is phenomenologi- 
cally parameterized in terms of the total impedance Z(u>) 
0, and one obtains a = Z '(0) / (h / '2e 2 ) . Such a trans- 
mission line model directly explains the doubling of the 
end exponent, as in the bulk case one has effectively two 
resistances in parallel as compared to the end case. Ap- 
plying this model to MWNTs, however, one finds rather 
small values for a, and, in addition, the observed volt- 
age dependence cannot be reproduced ||. On the other 
hand, for a LL, the elementary excitations are harmonic, 
I(lo) is known, and hence P(E) can also be computed 
exactly. The resulting exponents can be written in terms 
of the standard interaction parameter g, and the ratio 
a cnd/ aLL i s ff-dependent and always larger than two ||H • 
In MWNTs g is modified by inner-shell screening, and the 
relation between a LL and g is affected by the available 
subbands |10| . Estimating a LL and a^ d for the situa- 
tion of Ref. H , they are at least one order of magnitude 
smaller than observed. Hence neither conventional CB 
nor LL theory can explain these data. 

For a ballistic system, the full crossover in the CB be- 
tween a (single-channel) LL and a (many-channel) classi- 
cal resistor has been worked out in Ref. jll| . The results 
of this study imply that it is crucial to take disorder into 
account here. We thus have to (i) establish the harmonic 
nature of the electromagnetic modes, (ii) compute the 
low-frequency spectrum I(u>) and hence the exponent a, 
and (iii) compute the end exponent. Note that a pertur- 
bative treatment of the interactions is not sufficient, as 
the power law (Q) is inconsistent with conventional (ID 
or 2D) Altshuler-Aronov-Lee (AAL) predictions for the 
TDOS jl2]. However, 2D AAL logarithmic dependencies 
have been observed for MWNT bundles, which presum- 
ably are characterized by weaker interactions |l3[ . 

The main difference between MWNTs and SWNTs, 
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apart from the larger radii of MWNTs, R w 4 to 10 
nm, stems from the presence of inner shells. Although 
inter-shell tunneling is largely suppressed for a number of 
reasons Jl4| , and hence transport proceeds only through 
the outermost contacted shell ||, inner shells cause a 
screening of the electron-electron interaction. For a com- 
putation of the TDOS, the latter effect, as well as spin 
and Fermi (K) point degeneracy Q, can be included by 
a renormalization of the electron-electron interaction po- 
tential U(q) p0| , and hence we consider spinless electrons 
with only one K point in what follows. Since different 
shells always have incommensurate lattices, a quasiperi- 
odic ionic potential from inner shells acts on outermost- 
shell electrons. The effect of such a potential is very 
similar to a random potential described by a mean free 
path £ — vt fl5|j , which can be estimated from a disor- 
dered tight-binding approach using random on-site ener- 
gies with variance <j\. For nearest-neighbor hopping t$ 
and 2N + 1 subbands (see below), one finds [[To) 



t = {t Q /<T E ) 2 R/(2N + 1) w lOi? 



(4) 



Here u E ~ to/10 corresponds to the hopping strength 
between adjacent shells, and N ps 5. Therefore disor- 
der is present even in impurity- free MWNTs. In ad- 
dition, "true" disorder imposed by imperfections, sub- 
strate inhomogeneities, or defects can be important. Us- 
ing Matthiessen's rule, £~ 1 — > t~ x + £^ , where £q is the 
mean free path due to "true" disorder. Typical values 
of I w 5 to 100 nm were reported in Refs. Hence 
disorder is still "weak", k F £ ^> 1, and the possibility of a 
diffusive phase of interacting electrons must be discussed. 
Doped MWNTs are in between the ID and 2D limits and 
may be expected to show diffusive behavior over a wide 
energy range, where the Anderson localization expected 
for truly ID disordered systems has not yet set in. We 
mention in passing that the phase relaxation length is 
£$ ps 500 nm at 2 K, with the usual T dependence due 
to electron-electron interactions . 

The bandstructure of a clean nanotube is described 
by a Dirac "light cone", E{k) = v\k\, with k = (k,k±) 
and quantized transverse momentum, fcj_ = n/R, where 
n = —N,...,N and N — [k F R]. The Fermi velocity is 
d = 8x 10 5 m/sec, and the number 2N + 1 of ID sub- 
bands arising from periodic boundary conditions around 
the circumference is determined by the doping level fi via 
kp = n/v. Since MWNTs studied experimentally so far 
are characterized by rather large doping, \x ps 0.5 eV jL6|, 
a typical value is N ps 5. The nth subband is character- 
ized by Fermi velocity and momentum, 



ly/i - {n/k F Rf 
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The theoretical description for small u is quite intricate 
and will be given elsewhere, as the disorder causes loga- 
rithmic divergences of the self energy due to the vanishing 



DOS of the Dirac cone p7[ . Fortunately, these complica- 
tions are absent for the case studied below, k F R ^> 1. In 
the following, we focus on energy scales E < v /R, where 
it is sufficient to take a fixed N and thereby ignore van 
Hove singularities associated with the opening of new ID 
bands as the energy varies. 

Let us start with the bulk TDOS, where standard 
diagrammatic perturbation theory |T^ ] adopted to the 
MWNT geometry is quite illuminating. The disorder- 
averaged Greens function is 

G R/A {E,k) = [E + v-E^+isgniE)^}- 1 , (6) 

and hence the non-interacting DOS is 



k i n=—N 
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Prefactors in expressions like Eq. (R) are always chosen to 
ensure the correct ID limit for N = 0. In the 2D limit, 
N ^> 1, the above DOS is related to the conventional 
2D DOS, i^2D = kp/^irv, via vq — 2ttRv2B- The vertex 
renormalization due to one impurity line is described by 
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where = m/R is also quantized. In the low-energy 
long- wavelength limit, ojt, \q\£ -C 1, Eq. (0) implies 



£(u, q) = 1 + ilut — Dq 2 r — D±q ± r , 
with diffusion constants D — (t/ttuq) ^ n v n and 



(9) 
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One checks easily that D± — > D for large kpR, with D — > 
v 2 t/2. For N = 0, however, = and D = v 2 t, and 
therefore both ID and 2D limits are correctly reproduced. 
Summing the usual ladder series then directly implies a 
diffusion pole, [1 — <f)] _1 . 

To ensure the validity of the ladder approximation, 
however, the smallness of "crossed" diagrams first needs 
to be checked. Since here one is actually in be- 
tween the ID and 2D limits, the usual reasoning lead- 
ing to a suppression factor l/k F £ for crossed diagrams 
does not apply. Evaluation of the simplest diagram 
of the crossed type and comparing it to the corre- 
sponding ladder diagram gives a suppression factor / = 
( 1 /2)E„l/(^o«™) 2 . Since / = 1/2 for JV = 0, it is 
clearly not justified to neglect crossed diagrams for small 
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doping level /i where one approaches the ID limit. The 
resulting breakdown of the diffusive picture in this limit 
is of course expected from ID localization theory. For 
UfR 1, however, we obtain / ~ In(kpR) /{2ir 2 kpR} <C 
1, and are thus entitled to neglect crossed diagrams. In 
that case, D±_ ~ D = v 2 t/2. 

Here we consider an effectively short-ranged ID inter- 
action, U{q) w Uo, since the 1/r tail of the Coulomb 
interaction potential is often externally screened, e.g. by 
close-by metallic gates. Even when working on an in- 
sulating substrate, a situation characterized by a long- 
ranged interaction, the finite length of the MWNT will 
lead to a cutoff for the logarithmic singularity. Note that 
Uq represents a ID Fourier component and hence is ap- 
propriate for screening lengths A s large compared to R. 
In the 2D limit, R S> A s , this interaction constant is re- 
lated to the respective 2D parameter via Uq = U2i)/R- 
For the case of interest here, R <C X s , and the order-of- 
magnitude estimate Uq/2hv 1 obtains implying a 
dielectric constant e ss 10. Under an extended Hubbard 
model description, this amounts to U/t ~ 1 to 5 in stan- 
dard terminology. The first-order interaction correction 
to the TDOS can then be obtained following well-known 
arguments |lj;|l^|. Since the Hartree correction is sub- 
leading, we focus on the exchange correction, 
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Together with Eq. (||) , this result describes the complete 
crossover from the ballistic LL to the diffusive limit. For 
r — > oo, no energy dependence is found, in agreement 
with a LL = to this order in Uo- In the diffusive limit, 
however, Eq. ( JTo| ) can be simplified to 

^ = -^Re V (^E + Dn 2 /R 2 Y 1/2 . 

(11) 

The n = contribution yields the E" 1 / 2 ID AAL correc- 
tion, while the n ^ contributions give a 2D logarithmic 
correction |12jl9l, 



Sv(E)/v =aln(E/E ) 



(12) 



with Eq a ju and a = U R/(2ttD). In fact, Eq. (|T|) also 
holds in the quasi-ballistic regime, Et 1 |Tq ] . Com- 
paring the ID and 2D contributions in Eq. (|ll|), we see 
that only for sufficiently low energy scales, E < E* , the 
ID AAL law, 5v ~ E^ 1 / 2 , takes over, where we estimate 



E* 



v( 
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< 1 meV 
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for the crossover scale. Equation ([l3]) holds for weak-to- 
intermediate interactions, but may change for very strong 



interactions. For E > E*, Eq. (|12j) then gives the leading 
contribution, and the TDOS should be essentially equiv- 
alent to the one of a 2D disordered interacting metal. 

Since in MWNTs, interactions are of intermediate 
strength, the lowest-order result ( |l2| ) is not sufficient. 
Higher-order terms in Uq were recently discussed by 
Kamenev and Andreev |l9). In the framework of the 
dynamical Keldysh formalism for a 2D system with long- 
ranged Coulomb interactions, a derivation of the TDOS 
valid for arbitrary interaction strength has been pre- 
sented. Within the saddle-point approximation for the 
emerging nonlinear a model, the electromagnetic modes 
are indeed Gaussian, with spectral density [jnj 



7T ^ (Da 2 - iuj) 2 1 
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Adopting this result to MWNTs in the energy regime 
E > E*, the bulk exponent a = I(oj — > 0) follows as 



UqR ln(l + vqUq) 



2nD 



VqUq 



(15) 



For weak interactions, the second factor is close to unity, 
and thus a is just the prefactor in Eq. (|l2]). In effect, 
the AAL logarithmic correction therefore exponentiates 
into a power law. Using Uq/2ttv w 1 and N ~ 5, we 
then estimate a « Rjl. To compare with experiment, 
note that t w lOi? from Eq. (^), so long as only intrinsic 
disorder via the inner-shell potentials is present. Then 
one gets the same exponent for different tubes, a ~ 0.1. 
When "true" disorder is included, a should increase and 
depend on the particular MWNT. This appears to be in 
qualitative agreement with Ref. 0] , where values between 
a = 0.24 and 0.37 were observed for 11 different MWNTs. 

Next we turn to the boundary TDOS by considering a 
semi-infinite MWNT, x > 0. In the absence of disorder, 
the boundary Greens function is (5 = + ) 



G b (E,k x ;x,x') = 



dk 



4 sin(fca;) sin(fcx') 



2?r E + n - E(k) + iS sgn(£) 



which can be written as Gq(x — x') — G' Q l (x + x') with the 
homogeneous Greens function Gq. Therefore the bound- 
ary simply causes the "image term" — Gq(x + x'). To 
include disorder, it is convenient to solve the Dyson equa- 
tion in a real-space formalism. In the Born approxima- 
tion, the self energy is S(x) = {ttvqt)- 1 ^^) -G%(2x)]. 
If we retain only the Gq(0) term in the self energy, the 
disorder-averaged boundary Greens function is 



G b (E,kx;x,x') 



dk 



— 4sin(fcx) sm(kx')G(E, k) , 



2tt 



(16) 
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where we use Eq. (g). In the ID limit, kpR< 1, the im- 
age part — Gq(2x) in the self energy is crucial and mod- 
ifies the structure of G b . For the more interesting limit 
fcp R ^> 1, it suffices to analyze a semi-infinite plane, with 
the lowest-order correction from the image part 



SG\r,X) 



1 



dx"G% 



x G%(2x" + 2X,0)G% 



(17) 



where the integration extends over the half plane, X — 
(x + x')/2 is the distance to the boundary, and f— x — x! 
is the relative coordinate. For r = 0, this can be seen to 
vanish identically. Under the conditions that kpX ^> 1 



and kpr 1, Eq. (17) simplifies to 



SG b (r,X) 



1 



kwlv< 



2L> 



-G h {2X)G h (r) 



This gives subleading corrections to Eq. ( [Tq ) that are 
suppressed by a factor 1/kpX or 1/kpr, whichever is 
smaller. As these corrections die out on the microscopic 
lengthscale 1/kp, we take Eq. ( |l6| ) in what follows. 

The real-space diagrammatic perturbation theory can 
then be carried out explicitly. Dressing the interaction 
vertices with impurity lines within the ladder approxima- 
tion, the first-order correction to the TDOS at position 
x in the diffusive limit is 



Sv(E,x) 



J q± 



gE,q) 



with £(E,q) given by Eq. (|). For x > (D/E) 1 / 2 , we 
recover the bulk result (]l0|). However, when tunneling 
into the end of the MWNT, x < (D/E) 1 / 2 , the pref- 
actor in Eq. (|l^ ) gets doubled. This doubling can be 
rationalized using the quasi-classical picture of Ref . [M , 
where the TDOS at position x is related to the Fourier 
transform of the return probability P{x,t). The latter 
satisfies the diffusion equation with the boundary condi- 
tion that no current flows across the boundary, leading 
to P(x,t)/P(oc,t) = 1 + exp(-x 2 /Dt). Close to the 
boundary, the return probability is then doubled. 

Repeating the above arguments leading to the power 
law (Q) for the bulk TDOS then implies again a power 
law, now characterized by the boundary exponent a cn d = 
2a. The doubling of the exponent as one goes from the 
bulk to the boundary is in agreement with the experi- 
mental results of Ref. ||. It should be noted that our 
derivation holds only in the diffusive limit. Remarkably, 
this doubling - which is based on the properties of the 
boundary diffuson close to the edge - coincides with the 
prediction of the classical resistor model. 

Let us conclude by listing several open questions raised 
by our study. Future work needs to address the situation 
at small doping levels, the magnetic field dependence of 



the TDOS, and the conductivity. We hope that future 
theoretical and experimental work continues to reveal the 
subtle and beautiful interplay of disorder, dimensionality, 
and interactions presented by multi-wall nanotubes. 
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